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Abstract The phase space analysis of cosmological parameters 17^ and 7 ^ is 
given. Based on this, the well-known quintessence cosmology is studied with an 
exponential potential y(0) = Vbexp(—A0). Given observational data, the current 
state of universe could be pinpointed in the phase diagrams, thus making the 
diagrams more informative. The scaling solution of quintessence usually is not 
supposed to give the cosmic accelerating expansion, but we prove it could educe 
the transient acceleration. We also hnd that the differential equations of system 
used widely in study of scalar field are incomplete, and then a numerical method 
is used to figure out the range of application. 
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1 Introduction 

The late time cosmic acceleration has been supported by many independent cos¬ 
mological observations, including the type la supernovae (SNIa) [ 1 ], large scale 
structure [2], cosmic microwave background (CMB) anisotropy [3] etc. Dark en¬ 
ergy is believed to have unveiled the mystery of the cosmic acceleration, via the 
ratio of its pressure and energy density as well as the equation of state (EoS) 
parameter w < —1/3. Depending whether EoS is constant or time-dependent, 
the dark energy candidate can be divided into many subclasses. The cosmological 
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constant model (ylCDM) with w = — 1 is the most robnst model bnt with the 
hne-tuning problem [il[5] and coincidence problem [ 6 ]. This led to a widespread 
speculation that the EoS of vacuum energy is not a constant. The quintessence 
model [7118119] was naturally introduced to study the possible existence of one or 
more scalar fields. 

The well-known quintessence scenario is described by a canonical scalar held 4> 
minimally coupled to gravity. The quintessence slowly evolving down a potential 
V{(p) can lead to the cosmic acceleration, which is similar to slow roll inhation 
in the early universe. The tracker solution of quintessence makes the universe to 
be insensitive to the initial conditions, thus alleviating the coincidence problem 
[ 6 ] . The dynamics of quintessence will have different characteristics under different 
potentials, and they have been studied in detail nniiiiiiiKis]. The phase space 
analysis is a most common and effective method, and has been used to investigate 
dynamical behavior of cosmological model, such as quintessence [mmi], tachyons 
[musi , fc-essence m, quintom [T^, phantom [19], modihed gravity Enum, and 
so on. In these literatures, the authors analyzed the dynamical evolution of a 
universe by the global attractor. However, their analyses are presented using some 
intermediate variables, such as x — y, whose definitions have no explicit physical 
significance; as a consequence, it is hard to discern the evolution of the universe. 

A transformation form x — y to — 7 ^ is more explicit in physical signif¬ 
icance, which also be used widely in study of scalar field. And then we try to 
analyse the evolution of cosmology in — 7 ^ space by employing this trans¬ 
formation. This cosmological parameters’ space shows the evolutionary history of 
the quintessence cosmology in a phase diagram with more explicit physical sig¬ 
nificance. Moreover, given observational data, more information can be obtained 
from the diagram. However the appearance of singularity indicates that the trans¬ 
formation is incomplete. Therefore, it needs to determine the scope of application 
for this transformation. 

This paper is organized as follows. In Sec. [21 the quintessence cosmology and 
the exponential potential are reviewed, and both forms of differential equations 
describing autonomous system are listed as well. In Sec. |3l two analytical methods 
based on two sets of differential equations are compared. In Sec. [H the singularity 
problem is described and discussed in detail. Finally, in Sec. |5l discussion and 
conclusion will be given. 


2 The quintessence cosmology 

Under the assumption of a flat and homogeneous universe, the Friedmann-Robertson- 
Walker (FRW) metric in units of SttG = c = 1 could be written as 



( 1 ) 


We consider a quintessence universe in the presence of a background prefect fluid 
which is described by EoS Wh = pt/pb ranging from 0 (dust) to 1/3 (radiation), 
where pb, Pb respectively represents the pressure and energy density of background 
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Table 1 The properties of the critical points (This table is taken from Ref. m)- 


fluid. The Friedmann equations of quintessence cosmology are 


= ^{Pb +Pcj,), ( 2 ) 

H =-^ {pb + P,l> + Pb + P4,), (3) 


where H is the Hubble parameter, is the energy density of scalar held, a dot 
represents the derivative with respect to cosmic time t. The pressure and energy 
density of quintessence are given by 


and 


= Y - 

(4) 

p<i> = ^ + ^{4>)- 

(5) 


The corresponding continuity equation of scalar field is 


dV 

. + 3770 + — = 0 . 
0(j) 


The EoS parameter of quintessence is 





( 6 ) 

(7) 


For V ^ (j)^ the quintessence approaches to cosmological constant model be¬ 
cause Wtf, ~ —1, while for V ^ 0^ the quintessence describes a stiff fluid with 


~ 1 . 


In order to transform the cosmological equations into an autonomous dynam¬ 
ical system, the most frequently used method is to introduce some auxiliary vari¬ 
ables |12| 


0 

^/6H’ 


y = 



1 dV 

vl4’ 


X = 


3H2’ 


( 8 ) 
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So Eqs. dll), ([3]) and ([ 6 ]) can be rewritten as 


/ /3. 2 ,3 , 2 2 , N 

X = Y -3/ - 1 ) 

+ ^Wbx{l-x'^ (9) 

/ -i _i_ ^ /I _i_ 2 2 n 

y =-y ^ 

+ ^W6y(l - a;^ - y^), (10) 

A'=-y6 A^(r-l)a;, (11) 


where 



and the prime denotes the derivative with respect to In a. 

Introducing the scalar field energy density and the effective equation of 
state 7 ,/, = 1 + as the following expressions 


1^0 

70 


P 0 2,2 

^=x +y , 

Prfi +P4) 


72 


P 0 


V + 02/2 


20^2 

a ;2 + 2/2 ’ 


Eqs. (l^- (fTT1) could be modified as [ 22 ] 


i70 = 3(76 - 70)170(1 - J70), 

70 = (2 - 70)(-370 + A^370f20), 
A' = -Y'370f20A^(r- 1). 


(13) 

(14) 


(15) 

(16) 
(17) 


Here, 76 = 1 + W6 is the effective EoS of background prefect fluid. The above 
system (ii5|)-(ini) becomes an autonomous system while the tracker parameter F 
is a function of the roll parameter A, which is widely used to study scalar field 
But it should be noted that the above equations are tenable only 
in a; > 0. If a; < 0, should be — In general, people use these 

equations just on the case a; > 0. 

The value of F depends on the form of potential V{<p), and it could be a 
constant when the potential 17(0) takes a special form, such as F = 0, 0.5,1 for 
17(0) = 17o0,17o0^, Vo exp(—A0), respectively. Here, Vo is a constant. In this paper, 
the phase space analysis of cosmological parameters and 70 is given. Based on 
this, the quintessence cosmology is studied with an exponential potential 17( 0 ) = 
17oexp(—A0). Under this potential, the parameter A is a constant and the above 
mentioned 3-dimensional autonomous system becomes 2-dimensional as a result. 

With regard to this exponential potential, the dynamical evolution was ana¬ 
lyzed by employing dynamical system techniques [I01[I3]. The critical points of the 
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Fig. 1 Phase space of x — y with A = 3 and W(, = 0. Point B is the global attractor. The red 
line is the unit circle which stands for = 1 +2/^ = !)• 

system and their properties are displayed in Table [T] Here w^ff is the effective 
EoS parameter, which is defined as 



P<l> ^ Pb Pb 
P4> + pb Pb P<1) + Pb 


2 2 , / -I 2 2 \ 

= X -y + Wb{l - X - y )■ 


(18) 


If Weff < —1/3, the universe will undergo an accelerated phase of expansion. In 
addition, Eqs. ©-dm) were used by the authors to plot the evolution diagram in 
the x — y parameter space. However, we find that it is more physical if we use Eqs. 
(II3-(ini) to plot the evolution diagram in the 7 ,^ — parameter space. In the 
following section, a comparison between these two approaches will be given. 

3 The comparison of two analyses 

There are totally Hve critical points in Table [1] Point O stands for a matter 
dominated universe (I7m = 1), but it is a saddle point for all values of A. PointsA± 
correspond to the universes dominated by the scalar field kinetic energy with no 
acceleration, and they are not the stable points. Points B and C are the more 
interesting critical points. And then, we will discuss their characteristics in detail 
in phase diagrams. 

3.1 Solution B 

In the case of point B, the scalar field energy density is 0 < = 3(1+Wb)/A^ < 1, 

implying the universe evolves under both the matter and scalar field. This solution 
is the famous scaling solution where the effective EoS matches the background EoS 
{Weff = Wb). In Fig.[Tl point B is plotted in the x — y parameter space, with A = 3 
and Wb = 0 (dust). Based on the result, point B is shown to be a global attractor, 
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Fig. 2 Phase space of — 7 ^ with A = 3 and 205 = 0. The blue point N stands for the 
current state of universe, and the black line is its evolution trajectory with ln(a) from —7 to 
fO. The green horizontal line represents = 2/3. The region below it stands for accelerated 
expansion of universe. 



since the system eventually evolves to point B regardless of the initial conditions. 
In addition, solution A- (or A+ when A is negative) represents the early universe 
as the arrows originate from it. However, the physical significance oi x — y plane 
is not enough transparent to describe current universe. Conventionally, solution 
point B is not considered as the representative of current universe, because the EoS 
of this point, Wb, cannot lead to an accelerating expansion. In order to alleviate 
this problem, the double exponential potential has been considered in Refs. |26l 
|9] . But we will see that this solution can actually provide the possibility of cosmic 
acceleration. 

On the contrast, a phase diagram plotted in 17^ — space can provide more 
physical information. Taking Fig. [2] as an example, the phase space analysis is 
performed and point B is shown to be the global attractor still, consistent with 
the result of Fig. [1] More importantly, the evolution of the universe could be 
presented in the — 7 ,^ phase space due to the cosmological context of the phase 
diagram. According to the latest observation m, it favors the current state of 
universe with = 0.685,7,^0 = 0.05 , which is labeled as the blue point N (Fig. 
[2]). Accordingly, the evolution curve of the universe is obtained using and 7,^0 
as the initial conditions. The phase space analysis in — 7 ^ plane indicates that 
the universe would evolve along the black line before reaching the stable point B. 
As mentioned before, the cosmic acceleration requires w^ff < —1/3, and hence 
70 < 2/3. As shown in Fig. [2] Point N evolving to point B along with black line 
will across the green horizontal line, that means the cosmic accelerating expansion 
will slow down until decelerating forever (state B). That is to say this solution 
could describe the current accelerated expansion of the universe, merely it will 
slow down until decelerating. This phenomenon is so-called transient acceleration, 
which is also suggested by SNIa data [28], and investigated in several theoretical 
scenarios [2mi30ll3n[32l[3^ . 

Very interestingly, the solutions of non-accelerating expansion are still found 
useful in presenting the current acceleration of the universe. This is because such 
solutions describe the final state of the evolution, leaving the intermediate evolu¬ 
tionary state of the accelerating expansion rather untouched. Therefore, the stable 
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Fig. 3 Same as Fig. ^but for A=0.5. 
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Fig. 4 Same as Fig. [2]but for A=0.5. 

point B is still valuable, since to trace back current cosmic accelerating expansion 
is now possible. 

On the other hand, the past evolution of the universe should be checked. The 
evolution with respect to ln(a) from —7 to 10 corresponding to redshift 2 ; from 
1095 to about —1 is calculated, which means the history of the universe from early 
epoch to the future should be tested. Starting from point N, the solid black line 
is found to be discontinuous in the past (Fig. [2]), implying a singularity in the 
process of time reversal. But it does not mean that the solution is invalid, and its 
reason will be discussed in the following Sec. SI 


3.2 Solution C 

In Table [U point C stands for the universe completely dominated by the scalar 
field = 1). If only matter (wb = 0) is considered, the effective EoS parameter 
is We// = W 41 = A^/3 — 1. For < 2, the exponential potential scalar field 
could drive the cosmological acceleration. In Fig. [3l the phase space analysis is 
performed for A = 0.5 and w;, = 0 in re — y parameter space. Point C is found to 
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Fig. 5 The initial conditions 7^0 A. The shaded red regions show their valid values. The 
perpendicular line stands for 7^0 = 2/3. 


be an attractor. The universe evolves from the initial state A- or A+ to the final 
state C. 

Fig. H] displays the phase space of universe in — 7 ^ parameter space with 
A = 0.5. Point C is now found to be the global attractor. Because of its dynamic 
behavior, solution C is insensitive to the initial conditions so that it could alleviate 
the coincidence problem. The universe evolves from decelerating to accelerating 
expansion. The current state, point N with = 0.685 and 7 <^o = 0.05, evolves 
to the global attractor C. The EoS parameter changes slightly, and the universe 
is expected to accelerate forever, which is like the cosmological constant. In the 
process of time reversal, point N could evolve to the unstable point A with no 
singularity. 


4 Discussion of singularity 

So far, the advantages of analysis in — 7 ^ space have been demonstrated. 
Nonetheless, the singularity problem just emerged in Sl^j, — 7^4 space, which is 
abnormal because these two analyses describe the same system. The reason is 
that, in some cases, the Eqs. (HSI-dlll) are no longer valid in the evolution of the 
system. As mentioned previously, the Eqs. (HB-dni) are established in a: > 0. In 
Eig.m any point in a; > 0 evolving to the past will enter into a; < 0 region which 
is not defined in — 7 ^ space, and then there will be a singularity in Fig. [2l 
For Fig. [3l the majority of point in x > 0 evolving to the past will not enter into 
X < 0 region, and then it does not have singularity in Fig. 31 This indicates that 
the Eqs. (I15I) - (I17I) can’t perfectly describe the dynamical system and they have a 
range of applications. However they have been used effectively in many studies, 
therefore, it should be necessary to determine their scope of application. Through 
the comparison of these two analyses, it shows that the value of A is crucial to the 
scope. We use numerical method to determine the appropriate value of A. 

To solve the differential Eqs. (na-dizi) of the system, the initial conditions 
700 and A are required. The 7^0 ~ A plane is scattered sufficiently for f/^o = 0.685. 
With these scattered points, the differential Eqs. (HSll-dni) could be solved in 
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Fig. 6 The evolution of solution B considering a: < 0. 


principle. Whether the effective and 'y 4 ,{z) have singularities will be used to 

check the validity of the initial conditions 7<^o and A. The result is shown in Fig. 
[5] with a shaded red region. The perpendicular line is 7^0 = 2/3, and the region 
to the left could lead to an accelerating expansion. As what could be seen, if the 
Eqs. (ll5ll - (fT7ll to be used, the initial value of A should be small, because bigger 
values of A usually lead to singularities in the cosmic evolution back to the past. 
Form the other perspective, if one moment of the past is chosen as initial point, 
the cosmological parameters evolving to the present will not pass the region in 
Fig.O 

From this perspective, the evolution of solution B described by 17^ — 7 ^ ap¬ 
pearing singularity does not mean this solution is invalid. It just has no definition 
in the past, but it can describe the future. If we consider the situation of j; < 0, 
the evolution of solution B is presented in Fig. ( 6 ] 


5 Conclusions and discussions 

Till now, the cosmological parameters’ space is employed, and the quintessence 
cosmology with an exponential potential V{4>) = Vbexp(—A</) is employed to il¬ 
lustrate the perspective. Previously, the scalar fields are analyzed using dynamical 
system techniques in the x—y parameter space. Now an analysis using a more phys¬ 
ically significant parameter space is made possible. By comparing these two 

analyses, the advantages of this cosmological parameters’ space are demonstrated. 
The information of system solutions could be visually described by its coordinate 
in the — 7^4 phase space. Given observational data, more information could 
be obtained. The solution B is the so-called scaling solution. Conventionally, this 
solution is not supposed to give the cosmic accelerating expansion, but we prove 
it could give a transient acceleration and this phenomenon has been supported 
by some observation data. In other words, the stable attractor solutions of the 
system are not the description of current state of the universe but the final state 
of cosmic evolution, while the intermediate evolutionary state remained uncertain. 
Unfortunately, there has been a singularity in evolution of solution B. Its reason is 
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that the transformation from x — y to — 7 ^ is incomplete, and then a numerical 
method is used to calculate the range of application. 

There are scaling solutions when in which the transient accelerating 

universe is given. While < 37 ;,, the universe is expected to accelerate forever, 
which is like the cosmological constant. In either case, the quintessence cosmology 
could describe the cosmic accelerating expansion in the present, moreover, it could 
alleviate the coincidence problem for its dynamic behavior. This demonstrates that 
the quintessence scalar field is a excellent scenario. In addition, a previous strong 
observational constraint implies A > 9 [3l]. Therefore, the transient accelerating 
universe is more expected. 
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